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Abstract 

We discuss cosmological solutions of Hofava-Witten theory describing the strongly coupled 
heterotic string. At energies below the grand-unified scale, the effective theory is five- not four- 
dimensional, where the additional coordinate parameterizes an S 1 /Z2 orbifold. Furthermore, it 
admits no homogeneous solutions. Rather, the static vacuum state, appropriate for a reduction 
to four-dimensional N = 1 supersymmetric models, is a BPS domain wall pair. Relevant 
cosmological solutions are those associated with this BPS state. In particular, such solutions 
must be inhomogeneous, depending on the orbifold coordinate as well as on time. We present 
two examples of this new type of cosmological solution, obtained by separation of variables 
rather than by exchange of the time and radius coordinates of a brane solution, as in previous 
work. The first example represents the analog of a rolling radii solution with the radii specifying 
the geometry of the domain wall pair. This is generalized in the second example to include a 
nontrivial "Ramond-Ramond" scalar. 
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1 Introduction 



Hof ava and Witten have shown that the strongly coupled Eg x Eg heterotic string can be identified 
as the eleven-dimensional limit of M-theory compactified on an S 1 jZ% orbifold with a set of E% gauge 
fields at each ten-dimensional orbifold fixed plane |2|. Furthermore, Witten has demonstrated 
that there exists a consistent compactification of this M-theory limit on a "deformed" Calabi-Yau 
threefold, leading to a supersymmetric N = 1 theory in four dimensions ||- Matching at tree 
level to the phenomenological gravitational and grand-unified (GUT) couplings [||, ||, one finds 
the orbifold must be larger than the Calabi-Yau radius, by a factor of ten or so. Since the GUT 
scale (about 10 16 GeV) is set by the size of the Calabi-Yau threefold, this implies that at energies 
below the unification scale there is a regime where the universe appears five-dimensional. This 
five-dimensional regime represents a new setting for early universe cosmology, which has been 
traditionally studied in the framework of the four-dimensional effective action. 

In a previous paper §, the effective five-dimensional Hof ava- Witten theory was derived for 
the universal fields, which are independent of the particular form of the Calabi-Yau manifold. It 
has a number of interesting and unusual features. The theory lives in a five-dimensional space 
which is a product of a smooth four-dimensional manifold times the orbifold S 1 jZ?,- As a result, 
it splits into a bulk N = 1, d = 5 supersymmetric theory with the gravity supermultiplet and 
the universal hypermultiplet, and two four-dimensional "boundary" theories which reside on the 
two orbifold fixed hyperplanes. The additional fields of the boundary theories are N = 1, d = 4 
gauge multiplets and chiral multiplets. The reduction from eleven-dimensional supergravity to five 
dimensions requires the inclusion of non-zero values of the four-form field strength in the internal 
Calabi-Yau directions. This leads to a gauged version of five-dimensional supergravity with a 
potential term that had not previously been constructed. In addition, the theory has boundary 
potentials for the projection of the bulk scalar field onto the orbifold planes. 

These potentials lead to a particularly interesting effect: the boundary sources mean that the 
theory has no solutions homogeneous in the orbifold direction. In particular, flat space is not a 
solution. Instead, the "vacuum" which leads to a supersymmetric flat four-dimensional space is a 
three-brane domain wall solution. The three-brane couples to the bulk potential and is supported 
by the sources provided by the two boundary potentials. More precisely, it is a double domain 
wall solution with the two 3 + 1-dimensional worldvolumes each covering an orbifold plane and the 
orbifold itself as the transverse coordinate. It is BPS, preserving half of the d = 5 supersymmetries, 
and so is the appropriate background for a further reduction to four-dimensional N = 1 supergravity 
theories. In such a reduction, four-dimensional space-time becomes identified with the three-brane 
world volume. 

Cosmologies in effective Hofava- Witten theory should thus be both five-dimensional and inho- 
mogeneous in the extra dimension. What should realistic models look like? In the ideal case, one 
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would have a situation in which the internal six-dimensional Calabi-Yau space and the orbifold 
evolve in time for a short period and then settle down to their "phenomenological" values while the 
three non-compact dimensions continue to expand. Then, for late time, when all physical scales 
are much larger than the orbifold size, the theory is effectively four-dimensional and should, in the 
"static" limit, provide a realistic supergravity model of particle physics. As we have argued above, 
such realistic supergravity models originate from a reduction of the five-dimensional theory on its 
domain wall vacuum state. Hence, in the "static" limit at late time, realistic cosmological solutions 
should reduce to the domain wall or a perhaps a modification thereof that incorporates breaking 
of the remaining four-dimensional N = 1 supersymmetry. Consequently, one is forced to look for 
solutions which depend on the orbifold coordinate as well as on time. The main goal of this paper 
is to present simple examples of such cosmological solutions in five-dimensional heterotic M-theory 
to illustrate some of the characteristic cosmological features of the theory. 

In earlier work J?], ||, we showed how a general class of cosmological solutions, that is, time- 
dependent solutions of the equations of motion that are homogeneous and isotropic in our physical 
d = 3 subspace, can be obtained in both superstring theories and M-theory defined in spacetimes 
without boundary. Loosely speaking, we showed that a cosmological solution could be obtained from 
any p-brane or D-brane by inverting the roles of the time and "radial" spatial coordinate. This 
method will clearly continue to work in Hofava-Witten theory as long as one exchanges time with 
a radial coordinate not aligned in the orbifold direction. An example of this in eleven-dimensions, 



based on the solution of |J, has been given in [10]. It can, however, not been applied to the 
fundamental domain wall since its radial direction coincides with the orbifold coordinate. This 
coordinate is bounded and cannot be turned into time. Also, as argued above, exchanging radius 
and time in the domain wall solution would not be desirable since it should be viewed as the 
vacuum state and hence should not be modified in such a way. Instead, the domain wall itself 
should be made time dependent thereby leading to solutions that depend on both time and the 
orbifold coordinate. As a result, we have to deal with coupled partial differential equations, but, 
under certain constraints, these can by solved by separation of variables, though the equations 
remain non- linear. Essentially, we are allowing the moduli describing the geometry of the domain 
wall and the excitations of other five-dimensional fields, to become time-dependent. Technically, 
we will simply take the usual Ansatze for the five-dimensional fields, but now allow the functions 
to depend on both the time and radial coordinates. We will further demand that these functions 
each factor into a purely time dependent piece and a purely radial dependent piece. This is not, 
in general, sufficient to separate the equations of motion. However, we will show that subject to 
certain constraints separation of variables is achieved. We can solve these separated equations and 
find new, cosmologically relevant solutions. In this paper, we will restrict our attention to two 
examples representing cosmological extensions of the pure BPS three-brane. A more general class 
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of solutions will be presented elsewhere [11]. 

The first example is simply the domain wall itself with two of its three moduli made time- 
dependent. We show that separation of variables occurs in this case. It turns out that these moduli 



behave like "rolling radii" [12] which constitute fundamental cosmological solutions in weakly cou- 
pled string theory. Unlike those rolling radii which represent scale factors of homogeneous, isotropic 
spaces, here they measure the separation of the two walls of the three-brane and its worldvolume 
size (which, at the same time, is the size of "our" three-dimensional universe). All in all, we there- 
fore have a time-dependent domain wall pair with its shape staying rigid but its size and separation 
evolving like rolling radii. 

For the second example, we consider a similar setting as for the first but, in addition, we allow 
a nonvanishing "Ramond-Ramond" scalar. This terminology is perhaps a little misleading, but 
relates to the fact that the scalar would be a type II Ramond-Ramond field in the case where 
the orbifold was replaced by a circle. This makes connection with type II cosmologies with non- 
trivial Ramond-Ramond fields discussed in [0, |§. Separation of variables occurs for a specific 
time-independent form of this scalar. The orbifold-dependent part then coincides with the domain 
wall with, however, the addition of the Ramond-Ramond scalar. This non-vanishing value of the 
scalar breaks supersymmetry even in the static limit. We find that the time-dependent part of 
the equations fits into the general scheme of M-theory cosmological solutions with form fields as 
presented in ref. 0, Applying the results of these papers, the domain wall moduli are found 
to behave like rolling radii asymptotically for early and late times. The evolution rates in these 
asymptotic regions are different and the transitions between them can be attributed to the nontrivial 
Ramond-Ramond scalar. 

Let us now summarize our conventions. We use coordinates x a with indices a, /3, 7, • • ■ = 
0, • • • , 3, 11 to parameterize the five-dimensional space M5. Throughout this paper, when we refer to 
the orbifold, we will work in the "upstairs" picture with the orbifold S 1 /Z2 in the x n -direction. We 
choose the range x 11 6 [— np, irp] with the endpoints being identified. The Z2 orbifold symmetry acts 
as x 11 — > —x 11 . Then there exist two four-dimensional hyperplanes fixed under the Z2 symmetry 
which we denote by M4 , i = 1,2. Locally, they are specified by the conditions x 11 = 0, np. The 
indices p, u, p, ■ ■ ■ = 0, • • • ,3 are used for the four-dimensional space orthogonal to the orbifold. 
Fields will be required to have a definite behaviour under the Z2 orbifold symmetry, so that a 
general field $ is either even or odd, with ^(x 11 ) = db$(— x 11 ). 



2 The five-dimensional effective action 



The five-dimensional effective action for Hof ava-Witten theory, obtained from the eleven-dimensional 
theory by compactifying on a Calabi-Yau three-fold, was derived in || for the universal zero modes; 
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that is, the five-dimensional graviton supermultiplet and the breathing mode of the Calabi-Yau 
space, along with its superpartners. These last fields form a hypermultiplet in five dimensions. 
Furthermore, the theory contains four-dimensional N = 1 gauge multiplets and chiral gauge mat- 
ter fields on the orbifold planes. To keep the discussion as simple as possible, we do not consider 
the latter. This simple framework suffices for the cosmological solutions we will study in this paper. 



The general Lagrangian will be presented elsewhere [13, 14 1. 

In detail, we have the five-dimensional gravity supermultiplet with the metric g a p and an Abelian 
gauge field A a as the bosonic fields. The bosonic fields in the universal hypermultiplet are the real 
scalar field V (the dilaton, measuring the volume of the internal Calabi-Yau space), the three- form 
C a/ 3ry and the complex Ramond-Ramond scalar £. Note that the three-form C a ^ can be dualized 
to a scalar field a. Hence the hypermultiplet contains four real scalar fields. As explained in the 
introduction, all bulk fields should be even or odd under the Z2 orbifold symmetry. One finds 
that the fields g^ u , flli,n, Ali, <r must be even whereas g^n, A^, C must be odd. If one studies 
cosmological solutions of the theory these transformation properties are important as they restrict 
the set of allowed solutions to those with the correct Z2 symmetry. Now consider the boundary 
theories. In the five-dimensional space M5, the orbifold fixed planes constitute the four-dimensional 

(i) 

hypersurfaces M\ , i = 1,2. With the standard embedding in the reduction from eleven to five 
dimensions, there will be an Eq gauge field AjP and gauge matter fields on the orbifold plane m[ . 

For simplicity, we will set these gauge matter fields to zero in the following. On the orbifold plane 

(2) (2) 
M4 there is an E$ gauge field . 

The five-dimensional effective action of Hof ava-Witten theory is then given by 

£5 = Sbulk + Abound (2-1) 

where 

Sbulk = - A / ^{R + \t<#F# + -^=^ 5e A a T Pl T 5e + -Ld a Vd a V 
Zk 5 Jm 5 I 1 y2 *v 

+ ^72 \ d <*° - *(f 9 «£ " " ^ e(x n )A a ] [d a a - i(£ff*£ - - 2a e(x u )A a ] 
+ ^d a m+^La%} (2.2) 

Abound =— r / ,„ v-gv "0 — 2 , n ,V-gv a 
4 JmP 4 Jm 2) 



1 2 r 

V / IvtrF^F^ -atrF^F^} . (2.3) 

167ra GU T ^ Jm® >- M ^ / V ) 



where T a p = d a Ap — dpA a and the F$ are the field strengths of the boundary gauge fields, 
while i?W' M/ = 7je fiupcr FpJ . Furthermore, K5 and ogut are the five-dimensional Newton constant 
and the gauge coupling respectively. The constant ag in the above action can be computed for 



a given internal Calabi-Yau space. Explicit formulae are presented in ref. J^, §]. In the above 
expression, we have dropped higher-derivative terms. The sigma-model for the scalar fields is the 
well-known coset Mq = SU(2, 1)/SU{2) x U(l) of the universal hypermultiplet. The coupling of 
a to A a implies that a U(l) symmetry on M.q has been gauged. This gauging also induces the 
ao-dependent potential term in (|2.3| ). It has been demonstrated fi [l4| that the above action is 
indeed the bosonic part of a minimal N = 1 gauged supergravity theory in five dimensions coupled 
to chiral boundary theories. 

The most striking features of this action from the viewpoint of cosmology (and otherwise) are 
the bulk and boundary potentials for the dilaton V in Skulk an d Sbound- These potential terms are 
proportional to the parameter «o an d their origin is directly related to the nonzero internal four- 
form that had to be included in the dimensional reduction from eleven dimensions. The boundary 
potentials lead to sources in the Einstein equation and the equation of motion for V and a that are 
proportional to 8(x 11 ) or 5(x n — np). Hence, as long as V is finite (the internal Calabi-Yau space 
is compact) purely time-dependent solutions of the theory do not exist as they could never cancel 
these delta-function sources. One is therefore led to always consider dependence on time and the 
orbifold coordinate x 11 . The presence of a bulk potential proportional to V~ 2 seems to indicate 
that the dilaton has a runaway behaviour and the internal space decompactifies at late time. This 
picture, however, is too naive in that it ignores the boundary potentials and the Z2 symmetries of 
the fields. In fact, as we will show, the correct static domain wall vacuum of the theory depends on 
the orbifold direction in a way so as to exactly cancel these potentials. Consequently, it is important 
to note that for cosmological solutions based on the domain wall the time-dependent scale factors 
do not feel the potential terms. 



3 The domain- wall vacuum solution 

In this section, we would like to review the static "vacuum" solution of the five-dimensional Hofava- 
Witten theory, as given in || . As argued in the introduction, this solution is the basis for physically 
relevant cosmological solutions. It is clear from the five-dimensional action given in the previous 
section, that flat spacetime is not a solution of the equations of motion. It is precluded from 
being a solution by the potential terms, both in the bulk and on the boundaries. If not flat space, 
what is the natural vacuum solution? To answer this, notice that the theory ([O]) has all of the 
prerequisites necessary for a three-brane solution to exist. Generally, in order to have a [D— 2)-brane 
in a D-dimensional theory, one needs to have a (D — l)-form field or, equivalently, a cosmological 
constant. This is familiar from the eight-brane |15| in the massive type II A supergravity in ten 
dimensions JIB], and has been systematically studied for theories in arbitrary dimension obtained 



by generalized (Scherk-Schwarz) dimensional reduction [17]. In our case, this cosmological term is 
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provided by the bulk potential term in the action (^J) , precisely the term that disallowed flat space 
as a solution. From the viewpoint of the bulk theory, we could have multi three-brane solutions 
with an arbitrary number of parallel branes located at various places in the x 11 direction. However, 
elementary brane solutions have singularities at the location of the branes, needing to be supported 
by source terms. The natural candidates for those source terms, in our case, are the boundary 
actions. This restricts the possible solutions to those representing a pair of parallel three-branes 
corresponding to the orbifold planes. This pair of domain walls can be viewed as the "vacuum" of 
the five-dimensional theory, in the sense that it provides the appropriate background for a reduction 
to the d = 4, N = 1 effective theory. 

From the above discussion, it is clear that in order to find a three-brane solution, we should 
start with the Ansatz 

dsl = aiyfdx^dx^^ + biyfdy 2 (3.1) 
V = V{y) 

where a and b are functions of y = x 11 and all other fields vanish. The general solution for this 



Ansatz, satisfying the equations of motion derived from action (2.1), is given by 

a = clqH 1 / 2 



b = b H 2 H = ^-a \y\+h (3.2) 

V = b H 3 

where ao, &o an d /io are constants. We note that the boundary source terms have fixed the form 
of the harmonic function H in the above solution. Without specific information about the sources, 
the function H would generically be glued together from an arbitrary number of linear pieces with 
slopes ±\/2ao/3. The edges of each piece would then indicate the location of the source terms. The 
necessity of matching the boundary sources at y = and irp, however, has forced us to consider 
only two such linear pieces, namely y £ [0, irp] and y £ [— irp, 0]. These pieces are glued together 
at y = and irp (recall here that we have identified irp and —irp). To see this explicitly, let us 
consider one of the equations of motion; specifically, the equation derived from the variation of g^ v . 
For the Ansatz in (3.1), this is given by 



a 



' 2 a' b' 1 V' 2 a 2 b 2 _ V2a b 



a a z a b 12 V z 18 V z 3 V 

where the prime denotes differentiation with respect to y. The term involving the delta functions 
arises from the stress energy on the boundary planes. Inserting the solution ( |3,2| ) in this equation, 
we have 

d 2 H = ^a (5(y)-5(y-irp)) (3.4) 
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which shows that the solution represents two parallel three-branes located at the orbifold planes. 
Using the five-dimensional super symmetry transformations presented in ref. ||, one can check that 
this solution indeed preserves four of the eight super symmetries of the theory. 

Let us discuss the meaning of this solution. As is apparent from the Ansatz ([3.1|), it has 3 + 1- 
dimensional Poincare invariance and, as just stated, it preserves four supercharges. Therefore, a 
dimensional reduction to four dimensions on this solution leads to an N = 1 supergravity theory. In 
fact, this is just the "physical" four-dimensional effective theory of strongly coupled heterotic string 
theory which is the starting point of low energy particle phenomenology. This has been explicitly 



demonstrated in ref. [13, 14]. The two parallel three-branes of the solution, separated by the bulk, 
are oriented in the four uncompactified space-time dimensions, and carry the physical low-energy 
gauge and matter fields. Therefore, from the low-energy point of view where the orbifold is not 
resolved, the three-brane worldvolume is identified with four-dimensional space-time. In this sense 
the Universe lives on the worldvolume of a three-brane. It is the purpose of the following sections 
to put this picture into the context of cosmology; that is, to make it dynamical. Consequently, 
we are looking for time dependent solutions based on the static domain wall which we have just 
presented. 

4 The domain-wall cosmological solution 

In this section, we will present a cosmological solution related to the static domain wall vacuum 
of the previous section. As discussed in ref. |7|, H), a convenient way to find such a solution is to 
use Ansatz ( |3.1| ) where the y = x 11 coordinate in the functions a, b and V is replaced by the time 
coordinate r. However, in Hofava-Witten theory the boundary planes preclude this from being a 
solution of the equations of motion, since it does not admit homogeneous solutions. To see this 
explicitly, let us consider the goo equation of motion, where we replace a(y) — ► a(r),b(y) — > (3(r) 
and V — > 7(t). We find that 

a 2 a/3 If al 1 V2a 1 

cr a p 12 7 Z 18 7^ 3 P7 

where the dot denotes differentiation with respect to r. Again, the term containing the delta 
functions arises from the boundary planes. It is clear that, because of the y-dependence introduced 
by the delta functions, this equation has no globally defined solution. The structure of equation 



(4.1) suggests that a solution might be found if one were to let functions a, b and V depend 
on both r and y coordinates. This would be acceptable from the point of view of cosmology, 
since any such solution would be homogeneous and isotropic in the spatial coordinates x m where 
m, n, r, ■ ■ ■ = 1, 2, 3. In fact, the previous Ansatz was too homogeneous, being independent of the 
y coordinate as well. Instead, we are interested in solutions where the inhomogeneous vacuum 
domain wall evolves in time. 
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We now construct a cosmological solution where all functions depend on both r and y. We 
start with the Ansatz 



ds\ = -N(T,y) 2 dT 2 + a(T,y) 2 dx m dx n r ]mn + b(r,y) 2 dy 2 (4.2) 
V = V(r,y) 

Note that we have introduced a separate function N into the purely temporal part of the metric. 
This Ansatz leads to equations of motion that mix the r and y variables in a complicated non-linear 
way. In order to solve this system of equations, we will try to separate the two variables. That is, 
we let 

N(r,y) = n(r)a{y) 
a(r,y) = a(r)a(y) 

b(r,y)=P(r)b(y) (4.3) 
V(T,y)=j(T)V(y) 

There are two properties of this Ansatz that we wish to point out. The first is that for n = a = 



{3 = 7 = 1 it becomes identical to (3.1). Secondly, we note that n can be chosen to be any function 
by performing a redefinition of the r variable. That is, we can think of n as being subject to a 
gauge transformation. There is no a priori reason to believe that separation of variables will lead 



to a solution of the equations of motion derived from the action (2.1). However, as we now show, 
there is indeed such a solution. It is instructive to present one of the equations of motion. With 
the above Ansatz, the goo equation of motion is given by [] 

a 2 (a" a' 2 a' V IV' 2 a 2 b 2 (3 2 V2 b , c , . r/ ..p\ 
72 — + -2 - 777/ + 77777 + 77 772 77 ~ ^r a 0T-My) ~ % " *P))- 



b 2 \ a a 2 ab l2V 2 l8V 2 -f 2 3 V K xyJ vy ^7 




a/3 If 
a 2 ' a f3 12 7 



+ -75-7T77I (4-4) 



Note that if we set n = a = = j = l this equation becomes identical to (|3.3j). Similarly, if we 



set a = b = V = 1 and take the gauge n = 1 this equation becomes the same as (4T). As is, the 
above equation does not separate. However, the obstruction to a separation of variables is the two 
terms proportional to ao- Note that both of these terms would be strictly functions of y only if we 
demanded that ft oc 7. Without loss of generality, one can take 

= 1- (4.5) 

We will, henceforth, assume this is the case. Note that this result is already indicated by the 
structure of integration constants (moduli) in the static domain wall solution ( |3.2| ). With this 



1 From now on, we denote by a, 6, V the y-dependent part of the Ansatz (4.3) 



condition, the left hand side of equation ( |4.4j ) is purely y dependent, whereas the right hand side 
is purely r dependent. Both sides must now equal the same constant which, for simplicity, we take 
to be zero. The equation obtained by setting the left hand side to zero is identical to the pure y 
equation fl3.3|). The equation for the pure r dependent functions is 



^! + «£_ll = (46) 
a 2 a/3 12 7 2 



Hence, separation of variables can be achieved for the goo equation by demanding that (^J^) is 
true. What is more remarkable is that, subject to the constraint that = 7, all the equations of 
motion separate. The pure y equations are identical to those of the previous section and, hence, 
the domain wall solution fl3.2|) remains valid as the y-dependant part of the solution. 
The full set of r equations is found to be 

a 2 a 8 1 7 2 , ,, „, 

2 5_ 2 «A + |_|» + ^ +2 «| + ii! = (4.8) 

a an p p n or a p 4j z 



a an a 2 12 7 



^ + 3^ + ^-2!-^ = (4.10) 
7 07 P7 7^ 717 



In these equations we have displayed (5 and 7 independently, for reasons to become apparent shortly. 
Of course, one must solve these equations subject to the condition that = 7. As a first attempt 
to solve these equations, it is most convenient to choose a gauge for which 

n = const (4-11) 

so that r becomes proportional to the comoving time t, since dt = n(r)dr. In such a gauge, the 
equations simplify considerably and we obtain the solution 

a = A\t-t \ p 

(3 = 7 = J B|t-to| 9 (4.12) 

where 



and A, B and to are arbitrary constants. We have therefore found a cosmological solution, based 



on the separation Ansatz (4.3), with the y-dependent part being identical to the domain wall 
solution ( |3.2| ) and the scale factors a, (3, 7 evolving according to the power laws (4.12). This means 
that the shape of the domain wall pair stays rigid while its size and the separation between the 
walls evolve in time. Specifically, a measures the size of the spatial domain wall worldvolume (the 
size of the three-dimensional universe), while (5 specifies the separation of the two walls (the size 
of the orbifold). Due to the separation constraint 7 = j3, the time evolution of the Calabi-Yau 
volume, specified by 7, is always tracking the orbifold. From this point of view, we are allowing 



two of the three moduli in (3^2), namely ao and bo, to become time-dependent. Since these moduli 
multiply the harmonic function H, it is then easy to see why a solution by separation of variables 
was appropriate. 

To understand the structure of the above solution, it is useful to rewrite its time dependent 
part in a more systematic way using the formalism developed in ref. (?], J8|. First, let us define new 
functions a, (3 and 7 by 



a 



7 



J : n 



(4.14) 



and introduce the vector notation 



(a*) 



$ 

\l J 



d = (di) 



1 

w 



(4.15) 



Note that the vector d specifies the dimensions of the various subspaces, where the entry d\ = 3 
is the spatial worldvolume dimensions, d<i = 1 is the orbifold dimension and we insert for the 
dilaton. On the "moduli space" spanned by a we introduce the metric 



Gij 



2(di5ij - 

G n i = 

36 , 



didj) 



(4.16) 



which in our case explicitly reads 



Furthermore, we define E by 



G = -12 



( 1 \ \ 

3 
0-3 



E 



n 



(4.17) 



e d-a e 3a+/3 



n 



(4.18) 
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The equations of motion ( [4.7|) -( Oo| ) can then be rewritten as 



1 -jr % 
-Ea Get = , 
2 



(4.19) 



It is straightforward to show that if we choose a gauge n = const, these two equations exactly 
reproduce the solution given in ( |4.12| ) and (|4.13j ). The importance of this reformulation of the 
equations of motion lies, however, in the fact that we now get solutions more easily by exploiting 
the gauge choice for n. For example, let us now choose the gauge 



n 



„d-a 



(4.20) 



Note that in this gauge E = 1. The reader can verify that this gauge choice greatly simplifies 
solving the equations. The result is that 



n 



67 = Ct + k\ 
(3 = (6 ± 4 v / 3) Ct + k 2 



(4.21) 



where C, k\ and k 2 are arbitrary constants. Of course this solution is completely equivalent to the 
previous one, eq. ( 4.12] ), but written in a different gauge. We will exploit this gauge freedom to 
effect in the next section. 

To discuss cosmological properties we define the Hubble parameters 



H=±a 
dt 



(4.22) 



where t is the comoving time. From ( 4.12j ) and ( |4.14| ) we easily find 

/ p \ 



H 



t-t 



P 



(4.23) 



Note that the powers p satisfy the constraints 

p T Gp = d-p=l. 



(4.24) 



These relations are characteristic for rolling radii solutions [12] which are fundamental cosmological 



solutions of weakly coupled heterotic string theory. Comparison of the equations of motion ( 4.19 ) 
indeed shows that the scale factors a behave like rolling radii. The original rolling radii solutions 
describe freely evolving scale factors of a product of homogeneous, isotropic spaces. In our case, the 
scale factors also evolve freely (since the time-dependent part of the equations of motion, obtained 
after separating variables, does not contain a potential) but they describe the time evolution of the 
domain wall. This also proves our earlier claim that the potential terms in the five-dimensional 
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action ( |2.1| ) do not directly influence the time-dependence but are canceled by the static domain 
wall part of the solution. 

Let us now be more specific about the cosmological properties of our solution. First note from 



eq. ( 4.23 ) that there exist two different types of time ranges, namely t < to and t > to- In the 
first case, which we call the (— ) branch, the evolution starts at t — > — oo and runs into a future 
curvature singularity @, f| at t = to- In the second case, called the (+) branch, we start out in a 
past curvature singularity at t = to and evolve toward t — > oo. In summary, we therefore have the 
branches 



t E 



[—00, to] (—) branch 
[to, +00] (+) branch 



(4.25) 



For both of these branches we have two options for the powers p, defined in eq. ( [4.23| ), corresponding 
to the two different signs in eq. ( [4.13| ). Numerically, we find 



PT - 



+.06 



Pi - 



+.48 



\ 



V 



-.45 
-.08 



(4.26) 



+.81 

for the upper and lower sign in ( 4.13; ) respectively. We recall that the three entries in these 
vectors specify the evolution powers for the spatial worldvolume of the three-brane, the domain 
wall separation and the Calabi-Yau volume. The expansion of the domain wall worldvolume has 
so far been measured in terms of the five-dimensional Einstein frame metric . This is also 
what the above numbers p\ reflect. Alternatively, one could measure this expansion with the 
four-dimensional Einstein frame metric gfy so that the curvature scalar on the worldvolume is 
canonically normalized. From the relation 



9 (4) 



(fii.ii) 



1/2 (5) 



we find that this modifies p\ to 



Pi = Pi + 



P2 



(4.27) 



(4.28) 



In the following, we will discuss both frames. We recall that the separation condition (3 = 7 
implies that the internal Calabi-Yau space always tracks the orbifold. In the discussion we can, 
therefore, concentrate on the spatial worldvolume and the orbifold, corresponding to the first and 
second entries in ( 4.26| ). Let us first consider the (— ) branch. In this branch t £ [— 00, to] and, 
hence, t — to is always negative. It follows from eq. ( [4.23| ) that a subspace will expand if its p 
component is negative and contract if it is positive. For the first set of powers p^ in eq. (4.26) 
both the worldvolume and the orbifold contract in the five-dimensional Einstein frame. The same 
conclusion holds in the four-dimensional Einstein frame. For the second set, pi, in both frames 
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the worldvolume contracts while the orbifold expands. Furthermore, since the Hubble parameter 
of the orbifold increases in time the orbifold undergoes superinflation. 

Now we turn to the (+) branch. In this branch t G [to, oo] and, hence, t — to is always positive. 
Consequently, a subspace expands for a positive component of p and contracts otherwise. In 
addition, since the absolute values of all powers p are smaller than one, an expansion is always 
subluminal. For the vector p-j the worldvolume and the orbifold expand in both frames. On the 
other hand, the vector p\ describes an expanding worldvolume and a contracting orbifold in both 
frames. This last solution perhaps corresponds most closely to our notion of the early universe. 



5 Cosmological solutions with Ramond forms 

Thus far, we have looked for both static and cosmological solutions where the form fields £, j\& 
and a have been set to zero. As discussed in previous papers J?], |8), turning on one or several such 
fields can drastically alter the solutions and their cosmological properties. Hence, we would like to 
explore cosmological solutions with such non-trivial fields. For clarity, in this paper we will restrict 
the discussion to turning on the Ramond-Ramond scalar £ only, postponing the general discussion 



to another publication [11] 



The Ansatz we will use is the following. For the metric and dilaton field, we choose 

ds\ = -N(T,y) 2 dT 2 +a(T,y) 2 dx m dx n Vmn + b(T,y) 2 dy 2 (5.1) 
V = V(r,y). 

For the £ field, we assume that £ = £(r, y) and, hence, the field strength F a = d a ^ is given by 

F = Y(T,y), F 5 = X(r,y). (5.2) 

All other components of F a vanish. Note that since £ is complex, both X and Y are complex. Once 
again, we will solve the equations of motion by separation of variables. That is, we let 

N(r,y) = n(r)N(y) 
a(r,y) = a{r)a(y) 

b(r,y) = (3(r)b(y) (5.3) 
V(T,y) = 7 (r)V(y) 

and 

X(r,y) = X (r)X(y) (5.4) 
Y(r,y) =4>{r)Y{y) . (5.5) 
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Note that, in addition to the £ field, we have also allowed for the possibility that N(y) ^ a{y). 
Again, there is no a priori reason to believe that a solution can be found by separation of variables. 
However, as above, there is indeed such a solution, although the constraints required to separate 
variables are more subtle. It is instructive to present one of the equations of motion. With the 
above Ansatz, the goo equation of motion becomes f] 

b :(S(y)-S(y-, P ))^ 




3 V K w/ w ' ' 7 J 

p(& aP_^f\_^\Xl\xl_^\Yl\^l 
n 2 [a 2 a(3 12 7 2 J 3b 2 V 7 3n 2 V 7 1 ' ' 



Note that if we set X = Y = and N = a this equation becomes identical to ( |4.4| ). We now see 
that there are two different types of obstructions to the separation of variables. The first type, 
which we encountered in the previous section, is in the two terms proportional to ato- Clearly, we 
can separate variables only if we demand that 

P = 1 (5-7) 

as we did previously. However, for non- vanishing X and Y this is not sufficient. The problem, of 
course, comes from the last two terms in ( |5.6D . There are a number of options one could try in 
order to separate variables in these terms. It is important to note that X and Y are not completely 
independent, but are related to each other by the integrability condition d T X(r, y) = d y Y(r, y). We 
find that, because of this condition, it is impossible to obtain a solution by separation of variables 
that has both X(r,y) and Y(r, y) non- vanishing. Now X(r,y), but not Y(t, y), can be made to 
vanish by taking £ = £(t); that is, £ is a function of r only. However we can find no solution by 
separation of variables under this circumstance. Thus, we are finally led to the choice £ = £(y). In 
this case Y(t, y) = and we can, without lose of generality, choose 

X = 1 • (5.8) 



At this point, the only obstruction to separation of variables in equation ( |5.6| ) is the next to last 
term, N 2 \X\ 2 /3b 2 V^/. Setting 7 = const is too restrictive, so we must demand that 

X = ^coe*%) (5.9) 

where Co is a non-zero but otherwise arbitrary real constant and 9{y) is an, as yet, undetermined 
phase. Putting this condition into the £ equation of motion 

dy \~bv x ) =0 (5 - 10) 



2 In the following, N, a, b, V denote the y-dependent part of the Ansatz (5.3) 
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we find that is a constant Oq and a oc with arbitrary coefficient. Note that the last condition 
is consistent with the static vacuum solution ([^). Inserting this result into the </05 equation of 
motion 



a 
a 



~N 



(1 




I.YL 

6 V 



(5.11) 



we learn that N oc a with arbitrary coefficient. Henceforth, we choose N = a which is consistent 
with the static vacuum solution ( |3.2| ). Inserting all of these results, the goo equation of motion now 
becomes 



a 



i2 



+ 



a'b' 1 V 
~a~J + 12 V 2 



18 V 2 



-Y^o^(S(y) -d(y-TT P )) 



n 2 \a 2 a (3 12 7 2 



r 2 1 

i~ (5-12) 

3 7 



Note that the left hand side is of the same form as the static vacuum equation (3.3). The effect of 
turning on the £ background is to add a purely r dependent piece to the right hand side. Putting 
these results into the remaining four equations of motion, we find that they too separate, with the 
left hand side being purely y dependent and the right hand side purely r dependent. Again, we 
find that in these equations the left hand sides are identical to those in the static vacuum equations 
and the effect of turning on £ is to add extra r dependent terms to the right hand sides. In each 
equation, both sides must now equal the same constant which, for simplicity, we take to be zero. 
The y equations for a, b and V thus obtained by setting the left hand side to zero are identical to 
the static vacuum equations. Hence, we have shown that 



N 



b 
V 
X 



a H 



b n H z 



b H 



XqH 



1/2 
2 
3 
3 



H 



— a |y| + ho 



(5.13) 



where Xq = CQe te °a 1 b^ 2 is an arbitrary constant. 



The r equations obtained by setting the right hand side to zero are the following. 



a 2 a (3 1 7 2 
a 2 + a]3 ~ I272 



n 



3 /? 2 7 



(5.14) 



.at 



a n 



13 0h 



a 2 



dp I7 2 



2- - 2 + '- - '-- + + 2-'- + - c z 

a an p p n or a p 4 7^ 



n 



P 2 1 



a an a 2 1 7 2 Cn n 2 
1 1 — h — = 

a an a 2 12 7 2 3 /3 2 7 



(5.15) 
(5.16) 



15 



7 



dry $j 7 



777 



- + 3^ + ^ - ^ '- - 2c; 



7 «7 /?7 7 



717 



n 







(5.17) 



In these equations we have, once again, displayed (3 and 7 independently, although they should 
be solved subject to the condition (3 = 7. Note that the above equations are similar to the r 
equations in the previous section, but each now has an additional term proportional to Cg. These 
extra terms considerably complicate finding a solution of the r equations. Here, however, is where 
the formalism introduced in the previous section becomes important. Defining d, (5 and 7 as in 
( |4.14j ), and a,E and G as in Q4.15j ), ( f4.18| ) and Q4.17| ) respectively, the equations ( |5.14| )-( |5"?L7| ) can 
be written in the form 



1 



£d" Gd + ^ET 1 ^ = , 



where the potential £7 is defined as 



with 



£7 = 2c^ 5 



(5.18) 



(5.19) 



( 6 ^ 

q= 

We can now exploit the gauge freedom of n to simplify these equations. Choose the gauge 



(5.20) 



(5.21) 



where d is defined in ( [4,15| ), Then E becomes proportional to the potential U so that the potential 
terms in ( p,18| ) turn into constants. Thanks to this simplification, the equations of motion can be 
integrated which leads to the general solution 0, || 



a = c In |n — r I + w In 



ST 



+ k 



(5.22) 



K T~1 — T 

where t\ is an arbitrary parameter which we take, without loss of generality, to be positive and 



s = sign(< q,q>) 



(5.23) 



<q,q> 

The scalar product is defined as < q, q >= (f r G~ l q. The vectors w and k are integration constants 
subject to the constraints 



q ■ w 
w T Gw 
q- k 



1 



ln(cg| <q,q> |) 
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(5.24) 



This solution is quite general in that it describes an arbitrary number of scale factors with equations 
of motion given by ( 5.1£j ). Let us now specify to our example. For G and q as given in eq. (4.17) 
and (5.20) we find that 



< q, q >= 1 



(5.25) 



hence s = 1, and 



/ \ 

-2 



(5.26) 



Recall that we must, in addition, demand that (3 = 7. Note that the last two components of c 
are consistent with this equality. We can also solve the constraints ( 5.24| ) subject to the condition 
(3 = 7. The result is 



w 



( W 3 + I \ 

Qw 3 

V ^ ) 



( k 3 + § In eg \ 

k 3 



(5.27) 



where 



1 \/3 
6 12 



u>3 = ~- ± — ( 5 - 28 ) 

and /C3 is arbitrary. We conclude that in the gauge specified by (|5.21| ), the solution is given by 

1 



a= (w 3 + ^)ln ( — - — 
b \Ti — r 

(3 = — 2 In |ri — t| + 6t/;3 In 



+ A; 3 + - In c 



n - r 



(5.29) 



7 = -^ln|Ti - r| + w 3 ln ( — - — ) + k 3 

3 \n-Tj 

with W3 as above. As a consequence of s = 1, the range for r is restricted to 



< r < ti 



(5.30) 



in this solution. Let us now summarize our result. We have found a cosmological solution with 
a nontrivial Ramond-Ramond scalar £ starting with the separation Ansatz (^^). To achieve sep- 
aration of variables we had to demand that (3 = 7, as previously, and that the Ramond-Ramond 
scalar depends on the orbifold coordinate but not on time. Then the orbifold dependent part of 
the solution is given by eq. ( |5,13| ) and is identical to the static domain wall solution with the ad- 



dition of the Ramond-Ramond scalar. The time dependent part, in the gauge (5.21), is specified 
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by eq. ( 5.29 ). Furthermore, we have found that the time-dependent part of the equations of mo- 
tion can be cast in a form familiar from cosmological solutions studied previously |7|, ||. Those 
solutions describe the evolution for scale factors of homogeneous, isotropic subspaces in the pres- 
ence of antisymmetric tensor fields and are, therefore, natural generalizations of the rolling radii 
solutions. Each antisymmetric tensor field introduces an exponential type potential similar to the 



one in eq. ( 5.19| ). For the case with only one nontrivial form field, the general solution could be 



found and is given by eq. (5.22). We have, therefore, constructed a strong coupling version of these 
generalized rolling radii solutions with a one-form field strength, where the radii now specify the 
domain wall geometry rather that the size of maximally symmetric subspaces. We stress that the 
potential U in the time-dependent equations of motion does not originate from the potentials in 



the action (2A) but from the nontrivial Ramond-Ramond scalar. The potentials in the action are 
canceled by the static domain wall part of the solution, as in the previous example. 

From the similarity to the known generalized rolling radii solutions, we can also directly infer 
some of the basic cosmological properties of our solution, using the results of ref. |7|, H. We expect 
the integration constants to split into two disjunct sets which lead to solutions in the (— ) branch, 
comoving time range t £ [— oo, to], and the (+) branch, comoving time range t G [to, oo], respectively. 
The (— ) branch ends in a future curvature singularity and the (+) branch starts in a past curvature 
singularity. In both branches the solutions behave like rolling radii solutions asymptotically; that 
is, at t — > — oo,to in the (— ) branch and at t — > to,oo in the (+) branch. The two asymptotic 
regions in both branches have different expansion properties in general and the transition between 
them can be attributed to the nontrivial form field. 

Let us now analyze this in more detail for our solution, following the method presented in 
ref. |7|, Q. First we should express our solution in terms of the comoving time t by integrating 
dt = n(r)dr. The gauge parameter n(r) is explicitly given by 

n = e^'Vi - Tp^-Vr -1 (5-31) 

where 

x = d-w, A = d-c. (5.32) 

Given this expression, the integration cannot easily be performed in general except in the asymptotic 
regions r — > 0, T\. These regions will turn out to be precisely the asymptotic rolling-radii limits. 
Therefore, for our purpose, it suffices to concentrate on those regions. Eq. ( 5.31| ) shows that the 



resulting range for the comoving time depends on the magnitude of A and x (note that A is a 
fixed number, for a given model, whereas x depends on the integration constants). It turns out 
that for all values of the integration constants we have either x<Aorx>0>A. This splits the 
space of integration constant into two disjunct sets corresponding to the (— ) and the (+) branch 
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as explained before. More precisely, we have the mapping 

[—00, to] for x < A < , 



t E 



-) branch 

[to, +00] for x > > A , (+) branch 



(5.33) 



where to is a finite arbitrary time (which can be different for the two branches). We recall that 
the range of r is < r < t\. The above result can be easily read off from the expression (5.31) 
for the gauge parameter. Performing the integration in the asymptotic region we can express r in 
terms of the comoving time and find the Hubble parameters, defined by eq. ( |4.22| ), and the powers 
p. Generally, we have 



x— A 



at r ~ 
at t ~ n 



(5.34) 



Note that, from the mapping flpg ), the expression at r ~ describes the evolution powers at 
t — ► —00 in the (— ) branch and at t — to in the (+) branch; that is, the evolution powers in the 
early asymptotic region. Correspondingly, the expression for r ~ n applies to the late asymptotic 
regions; that is, to t — to in the (— ) branch and to t —* 00 in the (+) branch. As before, these 
powers satisfy the rolling radii constraints (4.24). 

Let us now insert the explicit expression for d, w and c, eqs. (P^) , (|5T27D and (pT2(j|) , that 
specify our example into those formulae. First, from eq. ( |5.32| ) we find that 



a/3 



(5.35) 



Note that the space of integration constants just consists of two points in our case, represented by 
the two signs in the expression for x above. Clearly, from the criterion fl5.33| ) the upper sign leads 
to a solution in the (+) branch and the lower sign to a solution in the (— ) branch. In each branch 
we therefore have a uniquely determined solution. Using eq. ( [5.34 ) we can calculate the asymptotic 
evolution powers in the (— ) branch 



P-,t- 



( +.06 \ 
+.81 

Correspondingly, for the (+) branch we have 

/ +.48 \ 
-.45 



P-,t-*to 



( +-48 \ 
-.45 
-.08 



(5.36) 



P+,t-**o 



V 



-.08 



P+,t- 



( +.06 \ 
+.81 

V + 13 / 



(5.37) 



Note that these vectors are in fact the same as in the (— ) branch, with the time order being 
reversed. This happens because they are three conditions on the powers p that hold in both 
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branches, namely the two rolling radii constraints ( [4.24 ) and the separation constraint (3 = 7 



eq. ( p.7|) , which implies that ps = 6p2- Since two of these conditions are linear and one is quadratic, 
we expect at most two different solutions for p. As in the previous solution, the time variation of the 
Calabi-Yau volume (third entry) is tracking the orbifold variation (second entry) as a consequence 
of the separation condition and, hence, needs not to be discussed separately. The first entry gives 
the evolution power for the spatial worldvolume in the five-dimensional Einstein frame. For a 
conversion to the four-dimensional Einstein frame one should again apply eq. ( [4.28 ). It is clear 



from the above numbers, however, that this conversion does not change the qualitative behaviour 
of the worldvolume evolution in any of the cases. Having said this, let us first discuss the (— ) 
branch. At t — > — 00 the powers are positive and, hence, the worldvolume and the orbifold are 
contracting. The solution then undergoes the transition induced by the Ramond-Ramond scalar. 
Then at t ~ to the worldvolume is still contracting while the orbifold has turned into super inflating 
expansion. In the (+) branch we start out with a subluminally expanding worldvolume and a 
contracting orbifold at t ~ to- After the transition both subspaces have turned into subluminal 
expansion. 

6 Conclusion 

In this paper we have presented the first examples of cosmological solutions in five-dimensional 
Hof ava-Witten theory. They are physically relevant in that they are related to the exact BPS three- 
brane pair in five dimensions, whose D = 4 worldvolume theory exhibits N = 1 supersymmetry. 



A wider class of such cosmological solutions can be obtained and will be presented elsewhere [11]. 
We expect solutions of this type to provide the fundamental scaffolding for theories of the early 
universe derived from Hofava-Witten theory, but they are clearly not sufficient as they stand. The 
most notable deficiency is the fact that they are vacuum solutions, devoid of any matter, radiation 
or potential stress-energy. Inclusion of such stress-energy is essential to understand the behaviour 
of early universe cosmology. A study of its effect on the cosmology of Hofava-Witten theory is 



presently underway [18]. 
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